Comprehensive molecular characterisation of the complete mitogenome of Ergasilus tumidus and phylogenetic relationships of Copepoda inferred from mitogenomes.

Although parasitic copepods of the genus Ergasilus von Nordmann, 1832 are globally distributed parasites of fish, their phylogenetic relationships with other Copepoda are not clear, and the characteristics of their mitochondrial genomes (mitogenomes) are not thoroughly understood. The objective of this study was to address these knowledge gaps by sequencing the complete mitogenome of Ergasilus tumidus Markevich, 1940. The complete mitogenome (GenBank Acc. No. OQ596537) was 14,431 bp long and it comprised 13 protein-coding genes (PCGs), 22 tRNAs, two tRNAs, and two control regions (CRs). Phylogenetic analyses, conducted using concatenated nucleotide and amino acid sequences of 13 protein-coding genes, produced two partially incongruent topologies. While the order Calanoida was consistently resolved as the sister lineage to the other three orders, topological instability was observed in the relationships of the orders Cyclopoida, Siphonostomatoida and Harpacticoida. Siphonostomatoida clustered with Cyclopoida in the nucleotide-based phylogeny, but with Harpacticoida in the amino acid-based phylogeny. The latter topology conforms to the widely accepted relationships, but we speculate that the former topology is more likely to be the correct one. Our study provides a complete mitogenome sequence of E. tumidus, which helps us better understand the molecular evolution of the genus Ergasilus. Additionally, we suggest a different perspective on the controversial phylogenetic relationships among Siphonostomatoida, Cyclopoida and Harpacticoida, diverging from previously accepted views.

Strong mitonuclear discordance in the phylogeny of Neodermata and evolutionary rates of Polyopisthocotylea.

The genomic evolution of Polyopisthocotylea remains poorly understood in comparison to the remaining three classes of Neodermata: Monopisthocotylea, Cestoda, and Trematoda. Moreover, the evolutionary sequence of major events in the phylogeny of Neodermata remains unresolved. Herein we sequenced the mitogenome and transcriptome of the polyopisthocotylean Diplorchis sp., and conducted comparative evolutionary analyses using nuclear (nDNA) and mitochondrial (mtDNA) genomic datasets of Neodermata. We found strong mitonuclear discordance in the phylogeny of Neodermata. Polyopisthocotylea exhibited striking mitonuclear discordance in relative evolutionary rates: the fastest-evolving mtDNA in Neodermata and a comparatively slowly-evolving nDNA genome. This was largely attributable to its very long stem branch in mtDNA topologies, not exhibited by the nDNA data. We found indications that the fast evolution of mitochondrial genomes of Polyopisthocotylea may be driven both by relaxed purifying selection pressures and elevated levels of directional selection. We identified mitochondria-associated genes encoded in the nuclear genome: they exhibited unique evolutionary rates, but not correlated with the evolutionary rate of mtDNA, and there is no evidence for compensatory evolution (they evolved slower than the rest of the genome). Finally, there appears to exist an exceptionally large (≈6.3 kb) nuclear mitochondrial DNA segment (numt) in the nuclear genome of newly sequenced Diplorchis sp. A 3'-end segment of the 16S rRNA gene encoded by the numt was expressed, suggesting that this gene acquired novel, regulatory functions after the transposition to the nuclear genome. In conclusion, Polyopisthocotylea appears to be the lineage with the fastest-evolving mtDNA sequences among all of Bilateria, but most of the substitutions were accumulated deep in the evolutionary history of this lineage. As the nuclear genome does not exhibit a similar pattern, the circumstances underpinning this evolutionary phenomenon remain a mystery.

Geography, phylogeny and host switch drive the coevolution of parasitic Gyrodactylus flatworms and their hosts.

BACKGROUND: Gyrodactylus is a lineage of monogenean flatworm ectoparasites exhibiting many features that make them a suitable model to study the host-parasite coevolutionary dynamics. Previous coevolutionary studies of this lineage mainly relied on low-power datasets (a small number of samples and a single molecular marker) and (now) outdated algorithms. METHODS: To investigate the coevolutionary relationship of gyrodactylids and their fish hosts in high resolution, we used complete mitogenomes (including two newly sequenced Gyrodactylus species), a large number of species in the single-gene dataset, and four different coevolutionary algorithms. RESULTS: The overall coevolutionary fit between the parasites and hosts was consistently significant. Multiple indicators confirmed that gyrodactylids are generally highly host-specific parasites, but several species could parasitize either multiple (more than 5) or phylogenetically distant fish hosts. The molecular dating results indicated that gyrodactylids tend to evolve towards high host specificity. Speciation by host switch was identified as a more important speciation mode than co-speciation. Assuming that the ancestral host belonged to Cypriniformes, we inferred four major host switch events to non-Cypriniformes hosts (mostly Salmoniformes), all of which occurred deep in the evolutionary history. Despite their relative rarity, these events had strong macroevolutionary consequences for gyrodactylid diversity. For example, in our dataset, 57.28% of all studied gyrodactylids parasitized only non-Cypriniformes hosts, which implies that the evolutionary history of more than half of all included lineages could be traced back to these major host switch events. The geographical co-occurrence of fishes and gyrodactylids determined the host use by these gyrodactylids, and geography accounted for most of the phylogenetic signal in host use. CONCLUSIONS: Our findings suggest that the coevolution of Gyrodactylus flatworms and their hosts is largely driven by geography, phylogeny, and host switches.

The energy metabolism of Balantidium polyvacuolum inhabiting the hindgut of Xenocypris davidi.

Anaerobic parasitic ciliates are a specialized group of ciliates that are adapted to anoxic and oxygen-depleted habitats. Among them, Balantidium polyvacuolum, which inhabits the hindgut of Xenocyprinae fishes, has received very limited scientific attention, so the molecular mechanism of its adaptation to the digestive tract microenvironment is still unclear. In this study, transmission electron microscopy (TEM) and single-cell transcriptome analysis were used to uncover the metabolism of B. polyvacuolum. Starch granules, endosymbiotic bacteria, and multiple specialized mitochondrion-related organelles (MROs) of various shapes were observed. The MROs may have completely lost the electron transport chain (ETC) complexes I, III, IV, and V and only retained succinate dehydrogenase subunit A (SDHA) of complex II. The tricarboxylic acid (TCA) cycle was also incomplete. It can be inferred that the hypoxic intestinal environment has led to the specialization of the mitochondria in B. polyvacuolum. Moreover, carbohydrate-active enzymes (CAZymes), including carbohydrate esterases, enzymes with a carbohydrate-binding module, glycoside hydrolases, and glycosyltransferases, were identified, which may constitute evidence that B. polyvacuolum is able to digest carbohydrates and starch. These findings can improve our knowledge of the energy metabolism and adaptive mechanisms of B. polyvacuolum.

Mitogenomic evolutionary rates in bilateria are influenced by parasitic lifestyle and locomotory capacity.

The evidence that parasitic animals exhibit elevated mitogenomic evolutionary rates is inconsistent and limited to Arthropoda. Similarly, the evidence that mitogenomic evolution is faster in species with low locomotory capacity is limited to a handful of animal lineages. We hypothesised that these two variables are associated and that locomotory capacity is a major underlying factor driving the elevated rates in parasites. Here, we study the evolutionary rates of mitogenomes of 10,906 bilaterian species classified according to their locomotory capacity and parasitic/free-living life history. In Bilateria, evolutionary rates were by far the highest in endoparasites, much lower in ectoparasites with reduced locomotory capacity and free-living lineages with low locomotory capacity, followed by parasitoids, ectoparasites with high locomotory capacity, and finally micropredatory and free-living lineages. The life history categorisation (parasitism) explained ≈45%, locomotory capacity categorisation explained ≈39%, and together they explained ≈56% of the total variability in evolutionary rates of mitochondrial protein-coding genes in Bilateria. Our findings suggest that these two variables play major roles in calibrating the mitogenomic molecular clock in bilaterian animals.

Mitochondrial metabolism of the facultative parasite Chilodonella uncinata (Alveolata, Ciliophora).

BACKGROUND: Chilodonella uncinata is an aerobic ciliate capable of switching between being free-living and parasitic on fish fins and gills, causing tissue damage and host mortality. It is widely used as a model organism for genetic studies, but its mitochondrial metabolism has never been studied. Therefore, we aimed to describe the morphological features and metabolic characteristics of its mitochondria. METHODS: Fluorescence staining and transmission electron microscopy (TEM) were used to observe the morphology of mitochondria. Single-cell transcriptome data of C. uncinata were annotated by the Clusters of Orthologous Genes (COG) database. Meanwhile, the metabolic pathways were constructed based on the transcriptomes. The phylogenetic analysis was also made based on the sequenced cytochrome c oxidase subunit 1 (COX1) gene. RESULTS: Mitochondria were stained red using Mito-tracker Red staining and were stained slightly blue by DAPI dye. The cristae and double membrane structures of the mitochondria were observed by TEM. Besides, many lipid droplets were evenly distributed around the macronucleus. A total of 2594 unigenes were assigned to 23 functional classifications of COG. Mitochondrial metabolic pathways were depicted. The mitochondria contained enzymes for the complete tricarboxylic acid (TCA) cycle, fatty acid metabolism, amino acid metabolism, and cytochrome-based electron transport chain (ETC), but only partial enzymes involved in the iron-sulfur clusters (ISCs). CONCLUSIONS: Our results showed that C. uncinata possess typical mitochondria. Stored lipid droplets inside mitochondria may be the energy storage of C. uncinata that helps its transmission from a free-living to a parasitic lifestyle. These findings also have improved our knowledge of the mitochondrial metabolism of C. uncinata and increased the volume of molecular data for future studies of this facultative parasite.

Inverted base composition skews and discontinuous mitochondrial genome architecture evolution in the Enoplea (Nematoda).

BACKGROUND: Within the class Enoplea, the earliest-branching lineages in the phylum Nematoda, the relatively highly conserved ancestral mitochondrial architecture of Trichinellida is in stark contrast to the rapidly evolving architecture of Dorylaimida and Mermithida. To better understand the evolution of mitogenomic architecture in this lineage, we sequenced the mitogenome of a fish parasite Pseudocapillaria tomentosa (Trichinellida: Capillariidae) and compared it to all available enoplean mitogenomes. RESULTS: P. tomentosa exhibited highly reduced noncoding regions (the largest was 98 bp), and a unique base composition among the Enoplea. We attributed the latter to the inverted GC skew (0.08) in comparison to the ancestral skew in Trichinellidae (-0.43 to -0.37). Capillariidae, Trichuridae and Longidoridae (Dorylaimida) generally exhibited low negative or low positive skews (-0.1 to 0.1), whereas Mermithidae exhibited fully inverted low skews (0 to 0.05). This is indicative of inversions in the strand replication order or otherwise disrupted replication mechanism in the lineages with reduced/inverted skews. Among the Trichinellida, Trichinellidae and Trichuridae have almost perfectly conserved architecture, whereas Capillariidae exhibit multiple rearrangements of tRNA genes. In contrast, Mermithidae (Mermithida) and Longidoridae (Dorylaimida) exhibit almost no similarity to the ancestral architecture. CONCLUSIONS: Longidoridae exhibited more rearranged mitogenomic architecture than the hypervariable Mermithidae. Similar to the Chromadorea, the evolution of mitochondrial architecture in enoplean nematodes exhibits a strong discontinuity: lineages possessing a mostly conserved architecture over tens of millions of years are interspersed with lineages exhibiting architectural hypervariability. As Longidoridae also have some of the smallest metazoan mitochondrial genomes, they contradict the prediction that compact mitogenomes should be structurally stable. Lineages exhibiting inverted skews appear to represent the intermediate phase between the Trichinellidae (ancestral) and fully derived skews in Chromadorean mitogenomes (GC skews = 0.18 to 0.64). Multiple lines of evidence (CAT-GTR analysis in our study, a majority of previous mitogenomic results, and skew disruption scenarios) support the Dorylaimia split into two sister-clades: Dorylaimida + Mermithida and Trichinellida. However, skew inversions produce strong base composition biases, which can hamper phylogenetic and other evolutionary studies, so enoplean mitogenomes have to be used with utmost care in evolutionary studies.

Evolutionary history of inversions in directional mutational pressures in crustacean mitochondrial genomes: Implications for evolutionary studies.

Inversions of the origin of replication (ORI) in mitochondrial genomes produce asymmetrical mutational pressures that can cause strong base composition skews. Due to skews often being overlooked, the total number of crustacean lineages that underwent ORI events remains unknown. We analysed skews, cumulative skew plots, conserved sequence motifs, and mitochondrial architecture of all 965 available crustacean mitogenomes (699 unique species). We found indications of an ORI in 159 (22.7%) species, and mapped these to 23 ORI events: 16 identified with confidence and 7 putative (13 newly proposed, and for 5 we improved the resolution). Two ORIs occurred at or above the order level: Isopoda and Copepoda. Shifts in skew plots are not a precise tool for identifying the replication mechanism. We discuss how ORIs can produce mutational bursts in mitogenomes and show how these can interfere with various types of evolutionary studies. Phylogenetic analyses were plagued by artefactual clustering, and ORI lineages exhibited longer branches, a higher number of synonymous substitutions, higher mutational saturation, and higher compositional heterogeneity. ORI events also affected codon usage and protein properties. We discuss how this may have caused erroneous interpretation of data in previous studies that did not account for skew patterns.

Slow crabs - fast genomes: Locomotory capacity predicts skew magnitude in crustacean mitogenomes.

Base composition skews (G-C/G+C) of mitochondrial genomes are believed to be primarily driven by mutational pressure, which is positively correlated with metabolic rate. In marine animals, metabolic rate is also positively correlated with locomotory capacity. Given the central role of mitochondria in energy metabolism, we hypothesised that selection for locomotory capacity should be positively correlated with the strength of purifying selection (dN/dS), and thus be negatively correlated with the skew magnitude. Therefore, these two models assume diametrically opposite associations between the metabolic rate and skew magnitude: positive correlation in the prevailing paradigm, and negative in our working hypothesis. We examined correlations between the skew magnitude, metabolic rate, locomotory capacity, and several other variables previously associated with mitochondrial evolution on 287 crustacean mitogenomes. Weakly locomotory taxa had higher skew magnitude and ω (dN/dS) values, but not the gene order rearrangement rate. Skew and ω magnitudes were correlated. Multilevel regression analyses indicated that three competing variables, body size, gene order rearrangement rate, and effective population size, had negligible impacts on the skew magnitude. In most crustacean lineages selection for locomotory capacity appears to be the primary factor determining the skew magnitude. Contrary to the prevailing paradigm, this implies that adaptive selection outweighs nonadaptive selection (mutation pressure) in crustaceans. However, we found indications that effective population size (nonadaptive factor) may outweigh the impact of locomotory capacity in sessile crustaceans (Thecostraca). In conclusion, skew magnitude is a product of the interplay between adaptive and nonadaptive factors, the balance of which varies among lineages.

Morphological and molecular characterization of a new ciliate Nyctotheroides grimi n. sp. (Armophorea, Clevelandellida) from Chinese frogs.

A new species of clevelandellid ciliate, Nyctotheroides grimi n. sp., is described from the frog Fejervarya limnocharis. Light and scanning electron microscopy were used for the morphological studies, and the DNA encoding the SSU rRNA gene (SSU rDNA) and the ITS1-5.8S subunit rRNA-ITS2 region (ITS) were sequenced for genetic comparisons and phylogenetic analysis. The main distinctive morphological feature is a knob-like projection in the left-posterior end; other differential characters are the cell size, the length of the oral groove and the shape of the infundibulum. Nyctotheroides grimi possess an apical suture line in the left and right side of the anterior end and in the left side of the caudal end. In the phylogenetic analyses, the new species engroups with other Nyctotheroides species forming a monophyletic group. The high similarity in the SSU rDNA and ITS sequences between Nyctotheroides species suggests a relative recent divergence. The genetic data and the different host range support the separation of Nyctotheroides and Nyctotherus; however the morphological criterion based on the presence (in Nyctotheroides)/absence (in Nyctothterus) of an apical kinetal suture line should be modified to consider the presence of kinetal suture lines in the apical and/or the caudal left side in Nyctotheroides.

Architectural instability, inverted skews and mitochondrial phylogenomics of Isopoda: outgroup choice affects the long-branch attraction artefacts.

The majority strand of mitochondrial genomes of crustaceans usually exhibits negative GC skews. Most isopods exhibit an inversed strand asymmetry, believed to be a consequence of an inversion of the replication origin (ROI). Recently, we proposed that an additional ROI event in the common ancestor of Cymothoidae and Corallanidae families resulted in a double-inverted skew (negative GC), and that taxa with homoplastic skews cluster together in phylogenetic analyses (long-branch attraction, LBA). Herein, we further explore these hypotheses, for which we sequenced the mitogenome of Asotana magnifica (Cymothoidae), and tested whether our conclusions were biased by poor taxon sampling and inclusion of outgroups. (1) The new mitogenome also exhibits a double-inverted skew, which supports the hypothesis of an additional ROI event in the common ancestor of Cymothoidae and Corallanidae families. (2) It exhibits a unique gene order, which corroborates that isopods possess exceptionally destabilized mitogenomic architecture. (3) Improved taxonomic sampling failed to resolve skew-driven phylogenetic artefacts. (4) The use of a single outgroup exacerbated the LBA, whereas both the use of a large number of outgroups and complete exclusion of outgroups ameliorated it.

Redescription of Opalina triangulata (Heterokonta, Opalinea) from Fejervarya limnocharis based on morphological and molecular data.

Opalinids are a large group of anaerobic protists, mainly inhabiting the cloacae of amphibians (frogs and toads). The classification of this group has not been fully resolved, because of a lack of molecular information. Here, we give a redescription of Opalina triangulata Metcalf, 1923, collected from the rectum of the frog Fejervarya limnocharis, based on detailed morphological and molecular data. Our phylogenetic analyses confirmed the monophyly of Opalinata. Within it, Opalinea were monophyletic with O. triangulata and O. undulata as well as two Protoopalina species grouping together. Karotomorpha and Proteromonas did not group together confirming the paraphyly of Proteromonadea. Meanwhile, the ITS2 secondary structural similarities as well as G-C content revealed greater similitudes between Opalina species and P. lacertae than with Blastocystis hominis, which is in accordance with their position as sister clades in the SSU rDNA-based phylogenies.

Mitochondrial Architecture Rearrangements Produce Asymmetrical Nonadaptive Mutational Pressures That Subvert the Phylogenetic Reconstruction in Isopoda.

The phylogeny of Isopoda, a speciose order of crustaceans, remains unresolved, with different data sets (morphological, nuclear, mitochondrial) often producing starkly incongruent phylogenetic hypotheses. We hypothesized that extreme diversity in their life histories might be causing compositional heterogeneity/heterotachy in their mitochondrial genomes, and compromising the phylogenetic reconstruction. We tested the effects of different data sets (mitochondrial, nuclear, nucleotides, amino acids, concatenated genes, individual genes, gene orders), phylogenetic algorithms (assuming data homogeneity, heterogeneity, and heterotachy), and partitioning; and found that almost all of them produced unique topologies. As we also found that mitogenomes of Asellota and two Cymothoida families (Cymothoidae and Corallanidae) possess inversed base (GC) skew patterns in comparison to other isopods, we concluded that inverted skews cause long-branch attraction phylogenetic artifacts between these taxa. These asymmetrical skews are most likely driven by multiple independent inversions of origin of replication (i.e., nonadaptive mutational pressures). Although the PhyloBayes CAT-GTR algorithm managed to attenuate some of these artifacts (and outperform partitioning), mitochondrial data have limited applicability for reconstructing the phylogeny of Isopoda. Regardless of this, our analyses allowed us to propose solutions to some unresolved phylogenetic debates, and support Asellota are the most likely candidate for the basal isopod branch. As our findings show that architectural rearrangements might produce major compositional biases even on relatively short evolutionary timescales, the implications are that proving the suitability of data via composition skew analyses should be a prerequisite for every study that aims to use mitochondrial data for phylogenetic reconstruction, even among closely related taxa.

Epidemiology and identification of two species of Chilodonella affecting farmed fishes in China.

The genus Chilodonella includes free-living ciliated protozoa as well as pathogenic species for freshwater fish, with Chilodonella hexasticha and Chilodonella piscicola being the most important ones. These parasites cause outbreaks with high mortalities among farmed freshwater fishes with great economic losses. There are few reports of these species in China, and their identification has been based mostly on their morphological characteristics. In the present work, the parasites causing five outbreaks occurring in China between 2014 and 2017 have been identified by morphological and genetic analysis. We provide the first records of Ctenopharingodon idella and Siniperca chuatsi as hosts of C. hexasticha, and of Procypris rabaudi and Schizothorax wangchiachii as hosts of C. piscicola. There are no differences in the gross pathological findings produced by C. hexasticha and C. piscicola, consisting in desquamation and necrosis of epithelial cells in the skin and gills and in severe fusion of gill lamellae. However, both species differ in their geographic distribution: C. piscicola was found in farms located at altitudes over 1500 m above sea level and with a water temperature ≤18 °C, while C. hexasticha was found in farms located at altitudes under 50 m above sea level and with a water temperature ≥21 °C. Present results confirm that C. hexasticha and C. piscicola are two different species that can be differenced by their morphology; however, their biological variability may lead to erroneous identifications and the diagnosis should be preferably based in genetic analysis including nuclear LSU rDNA and mitochondrial SSU rDNA sequences.

Dietary Bile Salt Types Influence the Composition of Biliary Bile Acids and Gut Microbiota in Grass Carp.

Lipid metabolism can influence host's health. There is increasing evidence for interplay between two key regulating factors in lipid metabolism: bile acids (BAs) and gut microbiota. However, very little is known about how types of different diet-supplemented bile salts (BS) influence this interaction in vivo. We sought to explore these relationships using grass carp (Ctenopharyngodon idellus), which often suffers functional disorder of liver and gallbladder. We studied fluctuations of BAs in the gall and changes of microbial communities in the gut in response to seven different diets: five different BS, chelating BS agent, and control. The BS comprised two primary BS [sodium taurochololate (TCAS) and sodium taurochenodeoxycholate (TCDCAS)], sodium tauroursodeoxycholate (TUDCAS), and two secondary BS [sodium taurodeoxycholate (TDCAS) and sodium taurolithocholate (TLCAS)]. Supplementation of primary BS caused a more significant fluctuation of biliary BAs than secondary BS, and TCAS caused a more prominent increase than TCDCAS and TUDCAS. For the gut microbiota, primary BS tended to increase their diversity and induce community succession, secondary BS resulted in a higher firmicutes/bacteroidetes ratio, while TUDCAS had no significant effects. Changes of the gut microbiota triggered by different types of BS caused alteration in BAs biotransformation. Two-obesity-associated families, Lachnospiraceae and Ruminococcaceae were positively correlated with biliary cholic acid (CA), taurochenodeoxycholic acid (TCDCA), and deoxycholic acid (DCA). As both primary and secondary BS resulted in increased synthesis of toxic secondary Bas by the gut microbiota, future studies should pay closer attention to gut microbiota when considering BA treatment.

The complete mitochondrial genome of Cymothoa indica has a highly rearranged gene order and clusters at the very base of the Isopoda clade.

As a result of great diversity in life histories and a large number of described species, taxonomic and phylogenetic uncertainty permeates the entire crustacean order of Isopoda. Large molecular datasets capable of providing sufficiently high phylogenetic resolution, such as mitochondrial genomes (mitogenomes), are needed to infer their evolutionary history with confidence, but isopod mitogenomes remain remarkably poorly represented in public databases. We sequenced the complete mitogenome of Cymothoa indica, a species belonging to a family from which no mitochondrial genome was sequenced yet, Cymothoidae. The mitogenome (circular, 14484 bp, A+T = 63.8%) is highly compact, appears to be missing two tRNA genes (trnI and trnE), and exhibits a unique gene order with a large number of rearrangements. High compactness and the existence of palindromes indicate that the mechanism behind these rearrangements might be associated with linearization events in its evolutionary history, similar to those proposed for isopods from the Armadillidium genus (Oniscidea). Isopods might present an important model system to study the proposed discontinuity in the dynamics of mitochondrial genomic architecture evolution. Phylogenetic analyses (Bayesian Inference and Maximum Likelihood) conducted using nucleotide sequences of all mitochondrial genes resolved Oniscidea and Cymothoida suborders as paraphyletic. Cymothoa indica was resolved as a sister group (basal) to all remaining isopods, which challenges the accepted isopod phylogeny, where Cymothoida are the most derived, and Phreatoicidea the most basal isopod group. There is growing evidence that Cymothoida suborder might be split into two evolutionary distant clades, with parasitic species being the most basal split in the Isopoda clade, but a much larger amount of molecular resources carrying a high phylogenetic resolution will be needed to infer the remarkably complex evolutionary history of this group of animals with confidence.

Altered gut microbiota associated with intestinal disease in grass carp (Ctenopharyngodon idellus).

Gut microbiota plays a crucial importance in their host. Disturbance of the microbial structure and function is known to be associated with inflammatory intestinal disorders. Enteritis is a significant cause of high mortality in fish species, including grass carp (Ctenopharyngodon idellus). Study regarding the association between microbial alternations and enteritis in grass carp is still absent. In this study, changes in the gut microbiota of grass carp suffering from enteritis were investigated using NGS-based 16S rRNA sequencing. Six healthy and ten abnormal fish (showing reddening anus, red odiferous fluid accumulating in the abdominal capacity, and flatulence and haemorrhage in the intestine) were collected from a fish farm in Huanggang Fisheries Institute (Hubei, China). Our results revealed that the diversity, structure, and function of gut microbiota were significantly different between diseased and healthy fish (P < 0.05). Particularly, members of the genera Dechloromonas, Methylocaldum, Planctomyces, Rhodobacter, Caulobacter, Flavobacterium, and Pseudomonas were significantly increased in diseased fish compared with that in healthy fish (P < 0.05). Predicted function indicated that microbiota significantly changed the specific metabolic pathways (related to amino acid metabolism, xenobiotics biodegradation and metabolism, and carbohydrate metabolism) in diseased fish (P < 0.05). Taken together, our findings point out the association between changes of the gut microbiota and enteritis in grass carp, which provide basic information useful for diagnoses, prevention, and treatment of intestinal diseases occurring in cultured fish.

The complete mitochondrial genome of parasitic nematode Camallanus cotti: extreme discontinuity in the rate of mitogenomic architecture evolution within the Chromadorea class.

BACKGROUND: Complete mitochondrial genomes are much better suited for the taxonomic identification and phylogenetic studies of nematodes than morphology or traditionally-used molecular markers, but they remain unavailable for the entire Camallanidae family (Chromadorea). As the only published mitogenome in the Camallanina suborder (Dracunculoidea superfamily) exhibited a unique gene order, the other objective of this research was to study the evolution of mitochondrial architecture in the Spirurida order. Thus, we sequenced the complete mitogenome of the Camallanus cotti fish parasite and conducted structural and phylogenomic comparative analyses with all available Spirurida mitogenomes. RESULTS: The mitogenome is exceptionally large (17,901 bp) among the Chromadorea and, with 46 (pseudo-) genes, exhibits a unique architecture among nematodes. Six protein-coding genes (PCGs) and six tRNAs are duplicated. An additional (seventh) tRNA (Trp) was probably duplicated by the remolding of tRNA-Ser2 (missing). Two pairs of these duplicated PCGs might be functional; three were incomplete and one contained stop codons. Apart from Ala and Asp, all other duplicated tRNAs are conserved and probably functional. Only 19 unique tRNAs were found. Phylogenomic analysis included Gnathostomatidae (Spirurina) in the Camallanina suborder. CONCLUSIONS: Within the Nematoda, comparable PCG duplications were observed only in the enoplean Mermithidae family, but those result from mitochondrial recombination, whereas characteristics of the studied mitogenome suggest that likely rearrangement mechanisms are either a series of duplications, transpositions and random loss events, or duplication, fragmentation and subsequent reassembly of the mitogenome. We put forward a hypothesis that the evolution of mitogenomic architecture is extremely discontinuous, and that once a long period of stasis in gene order and content has been punctuated by a rearrangement event, such a destabilised mitogenome is much more likely to undergo subsequent rearrangement events, resulting in an exponentially accelerated evolutionary rate of mitogenomic rearrangements. Implications of this model are particularly important for the application of gene order similarity as an additive source of phylogenetic information. Chromadorean nematodes, and particularly Camallanina clade (with C. cotti as an example of extremely accelerated rate of rearrangements), might be a good model to further study this discontinuity in the dynamics of mitogenomic evolution.

Succession and Fermentation Products of Grass Carp (Ctenopharyngodon idellus) Hindgut Microbiota in Response to an Extreme Dietary Shift.

Dietary intake affects the structure and function of microbes in host intestine. However, the succession of gut microbiota in response to changes in macronutrient levels during a long period of time remains insufficiently studied. Here, we determined the succession and metabolic products of intestinal microbiota in grass carp (Ctenopharyngodon idellus) undergoing an abrupt and extreme diet change, from fish meal to Sudan grass (Sorghum sudanense). Grass carp hindgut microbiota responded rapidly to the diet shift, reaching a new equilibrium approximately within 11 days. In comparison to animal-diet samples, Bacteroides, Lachnospiraceae and Erysipelotrichaceae increased significantly while Cetobacterium decreased significantly in plant-diet samples. Cetobacterium was negatively correlated with Bacteroides, Lachnospiraceae and Erysipelotrichaceae, while Bacteroides was positively correlated with Lachnospiraceae. Predicted glycoside hydrolase and polysaccharide lyase genes in Bacteroides and Lachnospiraceae from the Carbohydrate-Active enZymes (CAZy) database might be involved in degradation of the plant cell wall polysaccharides. However, none of these enzymes was detected in the grass carp genome searched against dbCAN database. Additionally, a significant decrease of short chain fatty acids levels in plant-based samples was observed. Generally, our results suggest a rapid adaption of grass carp intestinal microbiota to dietary shift, and that microbiota are likely to play an indispensable role in nutrient turnover and fermentation.